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1. Introduction 



The PACS-CS Collaboration has been pushing Np = 2 + 1 full QCD simulations toward the 
physical point with the nonperturbatively 0(a)-improved Wilson quark action and the Iwasaki 
gauge action [jlj |[ ||]. As presented in a separate report [Q] in detail, sizable amount of hadron spec- 
trum data has been collected at j8 = 1.9 on a 32 3 x 64 lattice down to the pion mass m n = 210MeV 
using the domain-decomposed Hybrid Monte Carlo (DDHMC) algorithm []|] on the PACS-CS 
computer. In this report we examine the chiral behavior of the pseudoscalar meson masses and 
decay constants in comparison with the prediction of chiral perturbation theory (ChPT). In partic- 
ular, we focus on the three points: (i) signals for chiral logarithms, (ii) determination of low energy 
constants in the chiral lagrangian, (iii) determination of the physical point with the ChPT fit. 



2. Simulation parameters 

We employ the f?(a)-improved Wilson quark action with a nonperturbative improvement co- 
efficient csw = 1-715[||] and the Iwasaki gauge action at /$ = 1.9 on a 32 3 x 64 lattice. In Table |l] 
we summarize the hopping parameters (k^, fc s ) and the results for the pion masses and the un- 
renormalized quark masses. The latter is defined through the axial vector Ward-Takahashi identity 
(AWI) by 

awi (V 4 A; mp (Qf (0)) nu 



where P is the pseudoscalar density and A™ p is the nonperturbatively 0(a) -improved axial vector 
current^]]. 

For later use we also define the renormalized quark mass and the pseudoscalar meson decay 
constant in the continuum MS scheme as follows: 



z A (i+b A mAm 



( m AWI \ C s he 1 

f PS = 2KuoZ A [l+b A ^— U^V— • < 2 - 3 ) 
y uq J Cp]j m PS 

Here C S A p are the amplitudes extracted from the correlation functions (A mip (f)/ 3 (0)) and (P(t)P(O)) 
with the exponentially smeared source and the local sink, while C l P is from (P(t)P(O)) with the local 
source and the local sink. The renormalization factors Z A j> and the improvement coefficients b A j> 
are evaluated perturbatively up to one-loop level[|| |J]with the tadpole improvement. 



3. Results 

3.1 Comparison with the previous CP-PACS/JLQCD results 

We first compare our results with the previous CP-PACS/JLQCD results, both of which are 
obtained with the same quark and gauge actions at j8 = 1 .9 but on different lattice sizes: 32 3 x 64 for 
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Kud 



am Tl 



am 



AWI 

ud 



am 



AWI 



MD time 



0.13700 
0.13727 
0.13754 
0.13770 
0.13781 



0.13640 
0.13640 
0.13640 
0.13640 
0.13640 



0.32196(62) 
0.26190(66) 
0.18998(56) 
0.13591(88) 
0.08989(291) 



0.02800(20) 
0.01895(13) 
0.01020(11) 
0.00521(9) 
0.00227(16) 



0.04295(30) 
0.04061(18) 
0.03876(18) 
0.03767(10) 
0.03716(20) 



2000 
2000 
2500 
2000 
350 



0.13754 0.13660 0.17934(78) 0.00908(7) 0.03257(17) 



900 



Table 1: Hopping parameters for the up-down and the strange quarks together with the results of am K , 
am ud l an d am t W1 - MD time is the number of trajectories multiplied by the trajectory length. 

the former and 20 3 x 40 for the latter. For the hopping parameters only the combination (k" U( j, fc s ) = 
(0.13700,0.13640) is in common, which is the heaviest case in the PACS-CS results whereas it is 
the lightest one in the CP-PACS/JLQCD results. 

In Fig. [j]we plot (am^) 2 / (am^ 1 ) and f K / ' f n as a function of am** 1 with K s fixed at 0.13640. 
The PACS-CS and the CP-PACS/JLQCD results are denoted by the black and the red symbols, 
respectively. The two sets of data together show a smooth behavior as a function of am^ 1 , and at 
fc ud = 0.13700 (am™ = 0.028) they show good consistency. 

An important observation is that while the CP-PACS/JLQCD results show an almost linear 
quark mass dependence both for {am n ) 2 / '(am™) and fx/fn, we find a clear curvature for the 
PACS-CS results, which is a characteristic feature of the ChPT prediction in the small quark mass 
region. This curvature drives the PACS-CS results for fa/ fit close to the experimental value toward 
the physical point. 
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Figure 1: Comparison of the PACS-CS (black) and the CP-PACS/JLQCD (red) results for (nm^/H^ 1 ) 
(left) and f K / 'fn (right). K s is fixed at 0.13640. 



3.2 Chiral fit formulae 

In ChPT, the one-loop expressions for the pseudoscalar meson masses and the decay constants 
are given by JlO|] 

m\ = 2mB |l+^-^7 ? + ^(16m(2L 8 -L 5 ) + 16(2m + m s )(2L 6 -L 4 ))|, (3.1) 
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m\ = (m + mjBo|l + ^ + ^(8(/n+m s )(2L 8 -L 5 ) + 16(2m+m s )(2L 6 -L4))| ; (3.2) 

U = fo 1 1 - 2\l* - Hk + j£ (8mL 5 + 8(2m + m s )L 4 ) | , (3.3) 
f 3 3 3 B 1 

/AT = /o 1 1 - ~ A Vn ~ ^ M# ~ ^Mt? + J2 + ' Ms ) Ls + 8 ( 2 ™ + m s) L *) | , (3.4) 

where m = (m u + m s )/2 and £4^6,8 are the l° w energy constants, and jUps is the chiral logarithm 
defined by 



1 mL 



with jit the renormalization scale. There are six unknown low energy constants So j/o> -£4,5,6,8 in the 
expressions above. The low energy constants are scale-dependent so as to cancel that of the chiral 



logarithm CiS). We determine these parameters by making a simultaneous fit for m 2 K , f n and 
fx- 

We also consider the contributions of the finite size effects based on ChPT At the one-loop 
level the finite size effect defined by Rx = (X(L) — X(°°))/X(°°) for X = mn-,mKifn,fK is given 
by [0]: 

R m n = \^ng\{^n)-^ n gl(\), (3.6) 

R mK = ^gi(^), (3.7) 
R fn = -tei{M-\$Kgi{kK), (3-8) 
Rh = -\S*8liM - \^m K ) - l^glilrj) (3.9) 



with 



m2ps Aps^mpsL, g l (x) = f i ^P-K l (^x), (3.10) 



(4*/,) 2 ' " ' °" ' „fi 

where T^i is the Bessel function of the second kind and m(n) denotes the multiplicities in the 

i] + n 2 z . 



expression of n = rr x + nt + n 2 . With the use of these formulae we estimate the possible finite size 



effects in our results. 



3.3 ChPT fits 



We apply the ChPT formulae (J3.l|)— (J3.4{) to our results at four points (j^ d , K&) = (0.13781, 
0.13640), (0.13770,0.13640), (0.13754,0.13640), (0.13754,0.13660). For these points, the p me- 
son mass satisfies the condition m p > 2m K . The measured AWI qurak masses are used for m and 



m s in eqs.(|j)-(0. The heaviest pion mass at (jc ud , K s ) = (0.13754,0.13640) is about 430MeV 
with the use of the cutoff determined below. The fit results are shown in Fig. ^, where the black 
solid lines are drawn with ?c s fixed at 0.13640 and the black dotted lines are for fc s = 0.13660. The 
red solid symbols represent the extrapolated values at the physical point whose determination is 
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L,(jU =m v ) 


PACS-CS 


PACS-CS with FSE 


exp. value[12|] 


MILC[13] 


U 


0.25(11) 


0.23(12) 


0.27 ± 0.8 


0.1(2)(2) 


L 5 


2.28(13) 


2.29(14) 


2.28 ±0.1 


2.0(3)(2) 


2L(, — L4 


0.16(4) 


0.16(4) 


0± 1.0 


0.5(1)(2) 


2U-L 5 


-0.59(5) 


-0.60(5) 


0.18 ±0.5 


-o.i(ixi) 


£ 2 /d.o.f. 


2.1(1.4) 


2.1(1.4) 







Table 2: Results for the low energy constants together with the phenomenological estimates and the MILC 
results. 



explained in Sec. 3.4. The heaviest point at (k^j, K s ) = (0.13754,0.13640) is not well described by 
ChPT both for (am„) 2 / (am^ 1 ) and fxj f%, and # 2 /d.o.f. is rather large (see Table [If). 

The results for the low energy constants are presented in Table where the phenomenological 



values with the experimental inputs[12] and the MILC results[13] are given for comparison. The 



renormalization scale is chosen to be = 0.547GeV. For L4 and L5 governing the behavior of 
f n , fx, our results show good agreement with both the phenomenological estimates and the MILC 
results. On the other hand, some discrepancies are observed between three results for 2L^ — L4 and 
2L% — L5 which enter into the ChPT formulae for m\ and m 2 K . 

In Fig. |2| we also draw the ChPT fit results incorporating the finite size effects. The blue 
solid lines are drawn for k s = 0.13640 and the blue dotted ones for k s = 0.13660. The fit curves 
with and without the finite size effects are almost degenerate for um^ 1 > 0.003, but we find a 
sizable difference at the physical point comparing the red open and solid symbols. This feature is 
understood by Fig. || where we plot the magnitude of Rx for X = m n ,mK,fK,fK with L = 2.8fm as 
a function of m K ( we note that R mps > and Rf ps < 0); the finite size effects are less than 2% for 
nips and fps at our simulation points. This is true even at the physical point except for f % which 
decreases by 4%. 



(am) /(am.) 




K =0.13640 

K =0.13660 

ChPT (K =0.13640) 

ChPT (K =0.13660) 

ChPT with FSE (K =0.13640) 

ChPT with FSE (K =0.13660) 



3.4 k 



f_/f 
K 71 




k =0.13640 

K =0.13660 

ChPT (K =0.13640) 

ChPT (K =0.13660) 

ChPT with FSE (K =0.13640) 

ChPT with FSE (K =0. 1 3660) 



0.005 



0.01 

AWI 

a m 



0.015 



0.02 



0.005 



0.01 

AWI 

a m 



0.015 



0.02 



Figure 2: Fit results for (ara,) 2 /(flm™) (left) and fx/fx (right). Red solid (open) symbols denote the 
extrapolated values at the physical point by the ChPT formulae without (with) the finite size effects. 



3.4 Physical point and light hadron spectrum 

In order to determine the up-down and the strange quark masses and the lattice cutoff we need 
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Figure 3: \Rx\ forX = in K ,mK,f7c,fK with L = 2.8fm as a function of m K . Solid vertical line denotes the 
physical point and the dotted ones are for our simulation points. 

three physical inputs. We try the following two cases: m^,m^,m^ and mx,mx,m<p. The choice of 
rriQ has theoretical and practical advantages: the Q. baryon is stable in the strong interactions and its 
mass, being composed of three strange quarks, is determined with good precision with small finite 
size effects. We also choose m$ for comparison. We employ the NLO ChPT formulae for the chiral 
extrapolations of m K , mx, f% and fx- A simple linear formula wihad = a + b- m^ 1 + c ■ mf^ 1 is used 
for the other hadron masses, employing data in the same range fc u d > 0.13754 as for pseudoscalar 
mesons. The results for the quark masses and the lattice cutoff are given by 

cT x =2.256(81)GeV, mjf = 2.37(1 l)MeV, mf* = 69.1(25)MeV, m^-input, (3.11) 
a" 1 =2.248(76)GeV, mjjf = 2.38(1 l)MeV, mf* = 69.4(25)MeV, m^-input, (3.12) 

where the errors are statistical. The two sets of results are consistent within the error. The quark 
masses are smaller than the recent estimates in the literature. We note, however, that we employed 
the perturbative renormalization factors to one-loop level which may contain a sizable uncertainty. 
A non-perutrbative calculation of the renormalization factor is in progress using the Schrodinger 
functional scheme. 

Using the physical quark masses and the cutoff determined above, we obtain predictions for 
the pseudoscalar meson decay constants at the physical point: 

f n = 144(6)MeV, f K = 175(6)MeV, f K /f n = 1.219(22), m n -input, (3.13) 
A = 143(6)MeV, = 175(5)MeV, fc/fn = 1-219(21), m -input, (3.14) 

to be compared with the experimantal values f n = 130.7 MeV, fx = 159.8 MeV, fx/ fn = 1.223. 
A 10% discrepancy in the magnitude of f n and fx might be due to use of one-loop perturbative 
since the ratio shows a good agreement. A non-nonperturbative calculation of Za and Z m = Za/Zp 
is also in progress. 

In Fig. ^| we compare the light hadron spectrum extrapolated to the physical point with the 
experiment. The results for the Q-input and the <p -input are consistent with each other, and both 
are in agreement with the experiment albeit errors are still not small for some of the hadrons. We 
find this to be encouraging. Further work is of course needed since cutoff errors of 0((aAQco) 2 ) 
are present in our results. 
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Figure 4: Light hadron spectrum extrapolated to the physical point with 12-input (red) and 0-input (blue). 
Horizontal bars denote the experimental values. 
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